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Abstract 

The mobility of graphene is very high because the quantum Hall effects 
can be observed even at room temperature. Graphene has the potential of 
the material for novel devices because of this high mobility. But the energy 
gap of graphene is zero, so graphene can not be applied to semiconductor 
devices such as transistors, LEDs, etc. In order to control the energy gaps, we 
propose atomic monolayers which consist of various materials besides carbon 
atoms. 

To examine the energy dispersions of atomic monolayers of various ma- 
terials, we calculated the electronic states of these atomic monolayers using 
density functional theory with structural optimizations. The quantum chem- 
ical calculation software Gaussian 03 was used under periodic boundary con- 
ditions. The calculation method is LSDA/6-311G(d,p), B3LYP/6-31G(d), 



* Corresponding author 
Email address: tatsuo@acp.metro-cit.ac.jp (Tatsuo Suzuki) 



or B3LYP/6-311G(d,p). 

The calculated materials are C (graphene), Si (silicene), Ge, SiC, GeC, 
GeSi, BN, BP, BAs, A1P, AlAs, GaP, and GaAs. These atomic monolayers 
can exist in the flat honeycomb shapes. The energy gaps of these atomic 
monolayers take various values. Ge is a semimetal; A1P, AlAs, GaP, and 
GaAs are indirect semiconductors; and others are direct semiconductors. 

We also calculated the change of energy dispersions accompanied by the 
substitution of the atoms. Our results suggest that the substitution of impu- 
rity atoms for monolayer materials can control the energy gaps of the atomic 
monolayers. 

We conclude that atomic monolayers of various materials have the poten- 
tial for novel devices. 



1. Introduction 

Graphene [l| is an atomic monolayer which consists of carbon atoms 
arranged in a honeycomb lattice. The mobility of graphene is very high 
because the quantum Hall effects can be observed even at room temperature 
0. Graphene has the potential of the material for novel devices because of 
this high mobility. But the energy gap of graphene is zero [If , so graphene can 
not be applied to semiconductor devices such as transistors, LEDs, etc. In 
order to control the energy gaps, we propose atomic monolayers which consist 
of various materials besides carbon atoms. Theoretical works show that the 
energy gaps of silicene which is the atomic monolayer of silicon atoms 0, 0, 111] 
and the atomic monolayer of germanium atoms [i| are zeros, but a pioneering 
experimental work about atomic monolayers of boron nitride formed on the 
metal surfaces \^ shows the finite energy gaps. So we hope that other atomic 
monolayers also have the finite energy gaps and have the potential for novel 
devices. To examine the energy dispersions of atomic monolayers of various 
materials, we calculated the electronic states of these atomic monolayers 
using density functional theory with structural optimizations. A honeycomb 
lattice was assumed as a crystal structure (Fig.l). The atoms were alternately 
placed like Fig. 1(b) in the case of a binary compound. 

2. Computational Method 

All calculations were performed using Gaussian 03 package [zj]. Gaussian 
is a very famous quantum chemical calculation software, and has been used 
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widely to calculate the electronic states of isolated molecules. Gaussian was 
improved in 2003, and Gaussian 03 has been able to calculate the electronic 
states using periodic boundary conditions. It can calculate the electronic 
states of not only three-dimensional crystals, but also two-dimensional or 
one-dimensional materials. 

We can specify the combination of method and basis set in Gaussian 
03 package. We want to calculate the electronic states of the Kohn-Sham 
equations of density functional theory, using the hybrid B3LYP exchange 
and correlation functionals combined with 6-311G(d,p) basis set (B3LYP/6- 
311G(d,p)). But unfortunately most calculations using B3LYP/6-311G(d,p) 
did not converge. So we used local spin density approximation (LSDA) com- 
bined with 6-311G(d,p) basis set (LSDA/6-311G(d,p)) or B3LYP combined 
with 6-31G(d) basis set (B3LYP/6-31G(d)) instead. The adopted combina- 
tions of method and basis set are noted in Table 1. 

We carried out the calculations of all electrons, so the contributions from 
core electrons are included. 

Furthermore the merits of our calculations are structural optimizations, 
and the translation vectors are also changed after every self-consistent cal- 
culation. 



3. Results and discussion 

3.1. A IV group element 

We show the energy dispersions of atomic monolayers which consist of a 
IV group element in the periodic table in Fig.2. We plot the lowest conduc- 
tion bands and the highest valence bands from other calculations by colored 
lines. 

Fig.2 (a) is a result of graphene which consists of carbon atoms. The 
conduction band and the valence band of graphene contacts at K points, 
which is famous as Dime points. The highest valence band of our calculation 
agrees very well with other calculations but the energy of the 



conduction band of our calculation at T point is rather higher than the results 
of other calculations. We cannot understand this reason at present. This is 
our problem to be solved. 

Fig.2 (b) is a result of silicene which consist of silicon atoms. The conduc- 
tion band and the valence band of silicene also contacts at K points as Dime 
points. Our result of silicene agrees very well with the solution of the local 
density approximation (LDA) with the projected augumented wave (PAW) 
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method [3j and the solution of LDA in a plane wave basis set with Vanderbilt 
ultrasoft pseudopotentials jij. The bond length of our calculation is 2.23 A, 
the bond length of LDA with the PAW method is 2.23 A [3], and the bond 
length of LDA with Vanderbilt ultrasoft pseudopotentials is 2.22 A Q. These 
bond lengths also agree very well. On the other hand, our result of silicene 
does not agree with the solutions of tight-binding model jij]. 

Fig. 2(c) is a result of an atomic monolayer which consists of germanium 
atoms. The conduction band and the valence band also contacts at K points, 
but the bottom of conduction band at T point is a little lower than the Fermi 
energy. So the atomic monolayer of germanium is a semimetal. This result 
agrees very well with the solution of LDA with the PAW method jjjj]. The 
bond length of our calculation is 2.31 A and the bond length of LDA with 
the PAW method is 2.32 A [3]. These bond lengths also agree very well. 

3.2. Binary compounds located in IV group elements 

We show the energy dispersions of atomic monolayers which consist of 
binary compounds located in IV group elements in the periodic table in Fig. 
3. 

The energy dispersions of SiC (Fig. 3 (a)) and GeC (Fig. 3(b)) resemble 
each other, and the conduction band and the valence band does not contact 
at K points. The energy gap of SiC is 3.526 eV and the energy gap of GeC 
is 3.160 eV. These values correspond to the wave lengths of ultraviolet rays. 

But the conduction band and the valence band of GeSi (Fig. 3(c)) almost 
contact at K points. This energy dispersion resembles the energy dispersion 
of silicene (Fig.2(b)). 

At present we cannot explain the physical reason why the conduction 
band and the valence band do not contact at K points for SiC and GeC 
while both bands contact for GeSi. 

3.3. Binary compounds located between III and V group elements 

We show the energy dispersions of atomic monolayers which consist of 
binary compounds located between III and V group elements in the periodic 
table in Fig. 4. We calculated A1N and GaN also, but the calculations of 
A1N and GaN did not converge. 

The energy gap of BN (Fig.4(a)) is very large (6.377 eV). There are 
experimental data of the valence band of atomic monolayers of BN formed 
on the metal surfaces 0]. Our result agrees well with this experimental band 
structures. 
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The energy dispersions of BP (Fig.4(b)) and BAs (Fig.4(c)) resemble 
each other. The energy gap of BP is 1.912 eV, which corresponds to the 
wave length of an orange light. The energy gap of BAs is 1.594 eV, which 
corresponds to the wave length of an infrared ray. 

The energy dispersions of A1P (Fig.4(d)), AlAs (Fig.4(e)), GaP (Fig.4(f)) 
and GaAs (Fig.4(g)) resemble one another. The bottoms of conduction bands 
are located at T point, so these are indirect semiconductors. 

We cannnot understand the physical factor which determines the transi- 
tion type (direct or indirect) at present. 

3. 4- Possibility of energy gap engineering 

We show that the change of energy dispersions accompanied by the sub- 
stitution of the atoms in Fig. 5. In the case of graphene (Fig.5(a)) and silicene 
(Fig. 5(e)), the conduction bands and the valence bands make Dirac points at 
K points, so the energy gaps are zeros. But when the 12.5% of carbon atoms 
are replaced with Si atoms (Fig.5(b)), the energy gap becomes open. This 
fact suggests that the substitution of impurity atoms for monolayer materials 
can control the energy gaps of the atomic monolayers. These atomic mono- 
layers are direct semiconductors, and the energy gap of SiC corresponds to 
the wave length of an ultraviolet ray. So these atomic monolayers can be 
applied to novel optical devices. 

We wanted to show the relation between the gap energy and the concen- 
tration of the substituted atoms. But by the difficulty of convergence, the cal- 
culation method is not unified. So we cannot estimate the gap-concentration 
relation quantitatively. 

3.5. Flatness 

In our all converged calculations, all atoms are located on the same plane 
which is formed by the two translation vectors. So the atomic monolayers of 
various materials can exist in the flat honeycomb shapes. 

4. Conclusion 

We calculated the electronic states of atomic monolayers which consist of 
various materials using density functional theory with structural optimiza- 
tions. The calculated materials are C (graphene), Si (silicene), Ge, SiC, GeC, 
GeSi, BN, BP, BAs, A1P, AlAs, GaP, and GaAs. These atomic monolayers 
can exist in the flat honeycomb shapes. The energy gaps of these atomic 
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monolayers take various values. Ge is a semimetal; A1P, AlAs, GaP, and 
GaAs are indirect semiconductors; and others are direct semiconductors. 

We also calculated the change of energy dispersions accompanied by the 
substitution of the atoms. Our results suggest that the substitution of impu- 
rity atoms for monolayer materials can control the energy gaps of the atomic 
monolayers. So we conclude that atomic monolayers of various materials 
have the potential for novel devices. 

Our important problem is to establish synthetic methods of these atomic 
monolayers. 
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Table 1: The results of calculations. 
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Figure 2: Energy dispersions of atomic monolayers which consist of a IV group element 
in the periodic table. The origin of energy is the Fermi energy. Black solid lines are 
our calculations, whose calcutation method is LSDA/6-311G(d,p). In Fig. 2(a), both 
red solid lines Q and red broken lines [9| are the local density approximation (LDA) 
with the GW approximation, green solid lines [HI are LDA with a full potential linear 
muffin orbital (FP-LMTO) method, and green broken lines ll| are LDA with a norm- 
conserving pseudopotential. In Fig. 2(b), red solid lines 0] are LDA with the projected 
augumented wave (PAW) method, red broken lines [H are LDA in a plane wave basis set 
with Vanderbilt ultrasoft pseudopotentials, green solid lines [B| are a first-nearest neighbor 
(INN) tight-binding model, and green broken lines la] are a second- nearest neighbor (2NN) 
tight-binding model. In Fig. 2(c), red solid lines |3| are LDA with the PAW method. 
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(a) SiC (b) GeC (c) GeSi 

Figure 3: Energy dispersions of atomic monolayers which consist of binary compounds 
located in IV group elements in the periodic table. The origin of energy is the Fermi 
energy. The calcutation method is B3LYP/6-31G(d). 
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Figure 4: Energy dispersions of atomic monolayers which consist of binary compounds 
located between III and V group elements in the periodic table. The origin of energy is 
the Fermi energy. The calcutation method of (a) is B3LYP/6-311G(d,p), and the method 
of others is B3LYP/6-31G(d). 
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(a) C (Graphene) (b) Si 125 C 875 ( c ) SiC 
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( d ) Si o.875 C o.i25 (e) Si (Silicene) 

Figure 5: (Top) The positions of atoms. Black balls are carbon atoms, and red balls are 
silicon atoms. Blue lines indicates an unit cell. (Bottom) Energy dispersions of atomic 
monolayers which consist of carbon atoms and silicon atoms. The origin of energy is 
the Fermi energy. The calcutation method of (a) and (e) is LSDA/6-311G(d,p), and the 
method of others is B3LYP/6-31G(d). 
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